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We discuss the path of a tracer particle as a microswimmer moves past on an infinite 
straight trajectory. If the tracer is sufficiently far from the path of the swimmer it moves 
in a closed loop. As the initial distance between the tracer and the path of the swimmer 
p decreases, the tracer is displaced a small distance backwards (relative to the direction 
of the swimmer velocity). For much smaller tracer-swimmer separations, however, the 
tracer displacement becomes positive and diverges as p — > 0. To quantify this behaviour 
we calculate the Darwin drift, the total volume swept out by a material sheet of tracers, 
initially perpendicular to the swimmer path, during the swimmer motion. We find that 
the drift can be written as the sum of a universal term which depends on the quadrupolar 
flow field of the swimmer, together with a non-universal contribution given by the sum of 
the volumes of the swimmer and its wake. The formula is compared to exact results for 
the squirmer model and to numerical calculations for a more realistic model swimmer. 
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1. Introduction 

Microswimmers, by virtue of their size, typically have Reynolds numbers <C 1. Hence 
they are well described by the Stokes equations. These are time independent and there- 
fore, as first recognised by Purcell (1977), any swimming stroke that leads to motion 
at Re=0 must be non-invariant under time reversal. Swimming strategies evolved by 
microorganisms to overcome this restriction include waving flagella and rotating helical 
filaments. 

Bacteria, algae and other microorganisms usually live in complex fluids, sharing their 
environment with a myriad of inert particles such as colloids and biofilaments. As they 
swim, the flow fields they produce stir the surrounding fluid, leading to the enhanced 



tracer diffusion observed in swimmer suspensions 


(Wu & Libchaber 


2000 Sokolov et al. 


2009 Leptos et al. 2009 Kurtuldu et al. 2011 


Miiio et al. 2011 


), and predicted by 


simulations and using theoretical arguments (Underhill et aL||2008 


Rushkin et aL||2010 


Ishikawa et al. 2010 Lin et al. 2011 Zaid et al. 


20111. Bacterial stirring is thought to 



be an important factor in controlling and enhancing nutrient uptake. Moreover, but still 
controversially, it has been suggested that stirring by small swimmers could provide a 
significant contribution to oceanic mixing ( [Katija fc Dabiri||2009| |Visser][2007l ) . 

f Email address for correspondence: mitya.pushkin@physics.ox.ac.uk 
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However, experimental trajectories of tracers in a suspension of swimmers often look 
loop-like (Leptos et al. 2009), with the net displacement between the initial and final 



locations considerably shorter than the characteristic trajectory size, Fig. [T] Indeed, it 
has been shown that tracers advected by a passing swimmer follow closed trajectories 
when the swimmer, modelled as a viscous dipole or quadrupole, travels along a straight, 
infinite path ( Dunkel et aL|20iO ). This loop-like nature of the swimmer paths was found, 
on the average, to survive the smearing effect of Brownian fluctuations. 

Clearly, closedncss of tracer trajectories severely limits the effectiveness of stirring, 
yet microswimmers enhance tracer diffusion. An explanation is called for, and this is an 
aim of this paper. We first, in Sec. 2, examine the tracer trajectories, identifying the 
properties of the flow field and the necessary assumptions and conditions that lead to 
loops. In particular we find that the trajectories are not closed for tracers near to the 
swimmer. Such particles are entrained by the swimmer, following it over a distance that 
diverges for a head-on collision. Loops are also destroyed by finite swimmer paths, and 
by the flow field appropriate to swimming in two dimensions. 

In the remainder of the paper we focus on investigating entrainment in detail. To do 
this we calculate the Darwin drift, the volume of fluid moved by a swimmer as it travels 
along an infinite straight path (Darwin 1953 Benjamin 1986 Eames et al. 1994). We 
find that although this quantity is infinite for a colloid pulled through a viscous fluid 
(Eames et al. 2003), it is finite for a swimmer. Moreover it can be written as the sum of 
a universal term, which depends only on the quadrupole flow produced by the swimmer, 
and a non-universal contribution equal to the volume of fluid in the swimmer wake. 

The analytic development leading to this result is described in Sec. 3 where we consider 
the displacement of a material sheet of tracers lying perpendicular to the swimmer path. 
In Sec. 3a we calculate the far field tracer displacement showing that it scales as p~ 3 for 
a dipolar swimmer, where p is the initial distance between the tracer and the swimmer 
path. We then, in Sec. 3b, obtain an expression for the Darwin drift produced by a 
passing swimmer. Section 4 compares the theoretical results to calculations for two model 
swimmers: In Sec. 4a we present analytical results for the squirmer model. Sec. 4b is 
devoted to a numerical calculation of the drift for a model bacterium which resembles R. 
sphaeroides. Algebraic and numerical details appear in Appendices. Finally our results 
are summarised, and other mechanisms for tracer diffusion discussed, in Sec. 5. 



2. Geometrical constraints to mixing by swimmers 

Consider a single swimmer moving steadily in a quiescent fluid at a speed V along 
the direction k. Let R(t) is the time-dependent position of the swimmer and u (r, k) the 
velocity induced at the position r relative to the swimmer. As the passive tracer moves 
with the local fluid velocity, its position satisfies 

^u(r T -R(i),k). (2.1) 

A self-propelled swimmer must be force-free. This physical requirement imposes a 
strong restriction on its velocity field. In particular, the Stokeslet term of the multipolar 
expansion must vanish and the velocity far from the swimmer is generally dominated by 
the force dipole flow: 



u(r,k) = -«(k-V)U s (r,k) ) (2.2) 

where Us is the Stokeslet, the gradient operator V acts on the first argument of Ug, 
and K is the dipole strength. The sign of k distinguishes the far-field of pushers (such 
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Figure 1 . Typical motion of a sheet of tracer particles as a swimmer moves in an infinite straight 
trajectory from left to right, along p = perpendicular to the sheet. The initial position of the 
tracer sheet is shown as a dotted line, and the envelope of the final tracer positions in blue. 
Examples of the loop like trajectories of the tracers are indicated as lighter, coloured lines. 
The starting and ending points of the tracer paths are marked by open circles and asterisks 
respectively and arrows indicate the direction of motion along the paths. For p > pc the net 
tracer displacement is negative and decays as p~ A as p — > oo for dipolar swimmers. For p < pc 
it is positive and diverges as p — > 0. However, the Darwin drift, the integrated volume swept out 
by the sheet, given by Eq (|3.9[), is finite. 



as E. coll), k > 0, from that of the pullers (such as C. reinhardtii) , k < 0. (Note that k 
needs to be multiplied by the prefactor 1/87777, where r\ is the kinematic fluid viscosity, 
in order to have units of a force dipole.) The Stokeslet can be expressed via the Oseen 
tensor J: 

U s (r,k)=k-J, J=^ + ^. (2.3) 

Here r = ||r|| and E is the identity matrix. Thus, the dipolar flow field u decays at large 
distances as r~ 2 . 

Note that the expression ( |2.2[ ) is considerably more restrictive than the general defi- 
nition of the force dipole (also commonly referred to as stresslet and Stokeslet doublet) 



(Chwang & Wu 19751. The strength of the general force dipole is a symmetric second 
order tensor and it may have six independent components. However, it is physically rea- 
sonable to expect that the lines of action of the resultant thrust and drag forces coincide 
and are aligned with the direction of swimming k. These conditions reduce the number 
of independent components from six to one. 

Under this assumption the quadrupolar terms of the multipole expansion for an ax- 
isymmetric swimmer can be written 

u« (r, k) = -1 (Q H (k • V) 2 + QxVi) U s (r, k) , (2.4) 



where Q\\ and Q\\ are the force quadrupole strengths along and perpendicular to the 

72^ 



swimming direction respectively and ~S/\ is the Laplace operator in the plane normal to 



the swimming direction. 
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In the particular case of Qn = Q± = Q, 



u 



Q 



(r,k)--|v 2 U 5 (r,k) = Q(k-V)V 



(2.5) 



A spherical body or a spherical droplet passively towed in viscous fluid has the quadrupole 
term of the same form ( Pozrikidis|1992 ). Remarkably, the flow field of a spherical particle 
moving in inviscid fluid also has the same form. Indeed, the expression to the right from 
the identity sign in (2.51 is but a potential dipole ( Pozrikidis]|1992 ) . 



It is essential for our further argument that the quadrupole term (2.4), similarly to the 



dipole term (2.2), can be expressed as a derivative along the swimming direction. Indeed, 
u Q (r,k) = (k- V)U Q (r,k), 

i\ Q±-Q\\ 



U Q (r,k) = QjV 



(k-V)Us (r,k) 



(2.6) 
(2.7) 



Hence a general flow field produced by a swimmer takes the form 

u (r, k) = (k • V) U (r, k) + O^). (2.8) 
We are now in a position to see that closedness of tracer trajectories is a universal 



feature rooted in the expression (2.8). Consider the velocity Uo(]"t — R-(i)) at the tracer 
location. Its Lagrangian derivative is 



dt 



= (V- V)U - 



dVT 

~d~t 



• V U , V = Vk. 



When the tracer velocity is negligible in comparison to the swimmer velocity, 
the Lagrangian derivative coincides with the Eulerian derivative: 



dr T 

dt I 



dUp 

dt 



V(k-V)U . 



Combining this expression with Eqns. (2.1) and (2.8), results in 

drx k dXJ, 

~dT ' 



a 



(2.9) 



(2.10) 



(2.11) 



V dt ' 

Therefore, the overall displacement of the tracer as the swimmer travels along the infinite 
path equals: 



drq 

~dt 



-dt 



V 



(U (+oo) - Uo(-oo)) 



0, 



(2.12) 



i.e. the tracer trajectory is closed. 

Besides translating, some swimmers, such as bacteria, often rotate around their bodies. 
But since typically the rotation axis coincides with the swimming direction, the induced 
azimuthal flow causes mixing only in the azimuthal direction. Our interest is in the cylin- 
drically symmetrical transport and hence we have disregarded the azimuthal component 



of the flow. The validity of the representation (2.2 ) was verified in recent experiments for 
E. coli in ( Drescher et aLpOlO ) and for C. reinhardtii in ( Guasto et a/.|[2010 ). 

The assumptions inherent in proving that the tracer loops are closed can be sum- 
marised as follows: 



1) The (induced) tracer velocity is much smaller than the swimmer velocity. 

2) The velocity field induced by the swimmer can be expressed in the form (2.8). 

3) The field U decays at infinity: U (r) — >• as r — > oo. 
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Table 1. Necessary conditions for closed tracer loops and corresponding mixing mechanisms. 



Closed loops when 



Examples of violation 



Mixing mechanism 



Induced tracer displacement 
<C Swimmer displacement 



Tracers close to swimmer 



Fluid entrainment 



u(r) = (k- V)Uo(r) 



Spinning bacterial cells 



Unimportant for 
axisymmetric swimmers 



Short-range interactions: 
||U (r)|| — > as r — > oo 



2D dipole 



Enhanced mixing in films 



Swimmers arrive from and 
leave to infinity 



Run-and-tumble of 

C. reinhardtii 



Random re-orientat ions (Lin 
et flt|2011[) 



4) Swimmers move in straights paths, arriving from and leaving to infinity. 



The effectiveness of biomixing ( Leptos et al.|2 009; Ku rtuldu et aL|2011 ) indicates that 
one or several of these assumptions are often violated in natural settings. The character 
of such violations may serve as basis for systematic analysis of different mechanisms of 
biomixing, and we list examples in table [l] In the remainder of this paper we provide 
a detailed analysis of situations where the first assumption is violated, and the tracer 
velocity is no longer small compared to that of the swimmer. This occurs for tracer 
particles that happen to lie in the vicinity of a swimmer path. Their motion can be 
understood as entrainment of the neighbouring fluid by a swimmer. Fuller comments on 
possible ways in which the other assumptions break down are given in the Discussion, 
Sec. 5. 



3. Fluid entrainment by swimmers 

The fluid flow near a swimmer is determined by the swimmer shape and gait. Con- 
sequently, important physical quantities, such as the dissipated energy or the ability to 
stir the surrounding fluid may strongly depend upon the particular swimmer features. 
However, some aspects of the hydrodynamics, for example the dipolar far flow field, are 
universal, i.e. generic across a wide range of swimmers. When considering the effect of 
entrainment on particle transport we shall, in particular, identify the universal aspects 
of the tracer motion. 

In order to analyse permanent displacement of fluid by a moving body it is convenient 
to envision an initially flat material sheet (ie a plane of tracer particles) of infinite exten- 
sion, and a swimmer approaching it normally, see Fig. [T] As the body moves, the sheet 
is deformed. Studies of the asymptotic shape of the material sheet and the fluid transfer 



fall under the rubric of the Darwin drift. The cases of towed bodies in inviscid (Darwin 
1953 Benjamin|1986 Eames et q/.|1994 | and viscous flows ( Eames et aZ.|2003 | have been 
analysed in detail. Here we consider the Darwin drift for force-free swimmers, far from, 
and then closer to, the swimmer trajectory. 
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3.1. Far field entrainment: 
The velocity field far from a swimmer is dominated by the dipolar and quadrupolar terms 



which, from Eqns. (2.2) and (2.4), are 

Q 



u (r, k) = -Kd z - - ^ Vi U S (r, k) + O 



(3.1) 



Here we use Cartesian coordinates with the z-axis directed along the swimming direction 
and hence (k • V) = d z . The multipole expansion ensures that the asymptotic shape of 
the material sheet is naturally universal at large distances from the swimmer path. Let 
p be the distance from this path to a tracer in the material sheet and A the resulting 
tracer displacement along the swimming direction, as the swimmer travels from and to 
infinity. Straightforward asymptotic calculations (see Appendix A) result in 



A 



-Ci 



V 2 p 3 



c, 



kQ± 1 



+ o( P - 5 ), 



for k^O. 



(3.2) 



Here C\ and are constants. Detailed calculation show that C\ — 7r/16. A typical 
shape of the deformed sheet is shown in Fig. [T] As C\ > 0, far from the swimmer the 
sheet is displaced oppositely to the swimming direction both for pushers and pullers. 



Pushers and pullers are distinguished by the sign of the second term in (3.2). Notably, 



the displacement A turns out independent of the quadrupole strength Qm . The calculation 
we have just described is valid in the region of small negative total tracer displacement. 
A quantitative comparison to numerical results will be given in Sec. 4b. 

3.2. Near field entrainment: 

Unlike the case of far distances, at close separations from the swimmer the shape of the 
deformed material sheet strongly depends on the swimmer size, shape, and stroke. It will, 
however, always remain integrable, and the volume of fluid entrained by the swimmer 
finite. We shall show, surprisingly, that the total volume of fluid displaced by the swimmer 
vd can be decomposed into a universal and a swimmer-dependent part having a lucid 
physical meaning. Before deriving this result, we pause to note that in his seminal paper 



Darwin ( 1953 ) put forward a remarkable proposal that for a body moving steadily in an 



inviscid fluid, vu equals the added mass of the body. He derived this result for a sphere 
and a cylinder of circular shape, finding in accordance with his hypothesis vd = CdVq, 
where Cd = 1/2 for a sphere and Cd = 1 for a cylinder. Later his arguments were refined 
( Eames et aL||199~4 and references therein) and a notion of the partial Darwin drift was 
introduced to account for the case of a body starting a finite distance from the sheet 
before crossing it and moving to infinity. More recently, a spherical viscous drop towed 



through a fluid of different viscosity was considered (Eames et al. 2003). It was found 



that in this case vd is divergent with a pre-factor proportional to the viscosity of the 



surrounding fluid. A viscosity-enhanced mechanism of ocean biomixing (Katija & Dabiri 



2009) was proposed based on this result. 



We follow the ideas of Darwin in our derivation of Vd for force-free swimmers. We 
set the z-axis along the swimming direction and introduce the cylindrically symmetric 
stream function in the system of reference co-moving with the swimmer tp(p,z). For a 



swimmer with the velocity field (3.1 ) we obtain (see Appendix B): 

'V 
2~ 



if) = 2irp 2 



KZ 



Q 



1-3 cos 2 9 



Q. 



l + 3cos 2 6»\ 1 



0{r~ 



(3.3) 




where r 2 = p 2 + z 2 and cos6> = z/r. 

Next, we consider a flow region around the swimmer which has a finite lateral extension 
—Lf z Li and is enclosed in a streamtube of a finite initial radius po : {p>z '■ 
ijj(p,z) tpQ, —Lf z ^ = ^i(po)-Li)}. The total entrained volume is obtained 

from the volume entrained by the flow within this region by letting po — > oo. By a simple 
geometrical argument (see Appendix B and Fig. [2]) , the latter equals 



v D = A(po) 



(3.4) 



Here is the volume of the flow region enclosed by a stagnation streamline around the 
swimmer (in the co-moving reference frame) and A depends only on the far-field of the 
swimmer velocity: 



2ttk 

A = — po sin ( 
tt(Q||+Q±) 



V 



cos ( 



+ 



AQ±-Q\\ 

V 



cos 



OiPo 1 )- 



(3.5) 



At this point, however, Eq. ( 3.5 1 is not well-defined, as the first term on the right hand side 



may diverge as po — > oo. But it will remain finite, in fact equal to zero, if the swimmer path 
is symmetric with respect to the initial plane of the material sheet, 8f = 7r — Hence, 
under this condition, the dipole flow component does not contribute to the permanent 
volume transfer across the sheet and, when the swimmer travels an infinite distance, 



A = 



M Q\\ + Q±) MQj--Q\\) 
~v + v 

The final expression for the Darwin drift reads 



V 



VD 



V 



(3.6) 



(3.7) 



Similarly to the tracer displacement A, the Darwin drift of a swimmer turns out inde- 



pendent of the quadrupole strength Qy. To appreciate the result (3.7), the distinction 



between the volume and the swimmer volume v s should be emphasized 



(3.8) 



where v wa f, e is the volume of the swimmer wake. 
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Eq. (3.7) provides a decomposition of vd into a universal and a swimmer-dependent 



part. For an arbitrary force- free swimmer the near-field details affect the entrained volume 



by virtue of only. In the next section we illustrate the utility of the expression (|3.7| 

mi 



for two model swimmers: the squirmer model and a numerical model of the bacterial cell. 
For the special case of Qn 



Q, the expression (3.7) reduces to 



vd = 



4ttQ 
V 



(3.9) 



This formula is, in fact, identical to the one derived by Taylor (1928) in the context 
of inviscid flows. But the underlying physical reasons are, of course, different. In the 



inviscid case the potential dipole (coinciding with the viscous force quadrupole (Chwang 
&: Wu|1975 )) is the only dominant fluid velocity component present that may contribute 
to A. In the (highly) viscous case, however, the Stokeslet and stresslet components do 
not contribute because of the force-free nature of the swimmer motion and the symmetry 
conditions on the swimmer trajectory respectively. We recently found that the expression 



(3.9) was also stated in the context of swimmers in Leshansky & Pismen (2010) but the 



questions of universality and the role of the wake were not discussed. 

We have assumed so far that microswimmers move in a roughly steady fashion and 
produce stationary flows. But some microswimmers, such as C. rheinhardtii, produce 
oscillatory flows by beating their fiagella asymmetrically during the power and recovery 
strokes. Nevertheless, all conclusions obtained in this chapter remain valid if we un- 
derstand the physical quantities such as the fluid velocity, the streamfunction and the 
multipole strengths in the beat cycle-averaged sense. 



4. Two model swimmers 

4.1. Entrainment by a spherical squirmer 
The spherical squirmer model has been the theorists workhorse for studying swimming 



at low Re numbers for more than half a century (Lighthill 1952 Blake 1971 Ishikawa] 
et al. 2006). It assumes a spherical swimmer of radius a with a prescribed tangential 
velocity ug on the swimmer's surface in the co-moving reference frame (9 is the polar 
angle measured from the swimming direction). The induced flow can be found analytically 
as a superposition of spherical harmonics. In the simplest nontrivial case when just the 
first two harmonics are nonzero: 

ug = Bi sin 9 + B 2 sin 9 cos 9, 



where B\ and B2 are constants. For a neutrally buoyant, force-free swimmer B\ sets 
the swimming speed V 
streamfunction equals 



the swimming speed V — \B\, and B2 sets the swimmer dipole strength. The squirmer 



irp 2 V 



1 



3/3 z 
~2r 



(*)' 



(4.1) 



where /3 = B\jBi- Comparison of the latter expression with (3.3) shows that squirmer 



dipole strength k = ^-Va 2 and the quadrupole strength Q± — Q\\ = Q = ^a 3 . Clearly, 
squirmers having /? > are pushers while those having < are pullers. For strong 
enough dipole strength, \/3\ > 1, squirmers develop a bubble- like wake. The dependence 
of the squirmer wake volume on /3, v wa k e /vo = /(/3), can be calculated analytically, see 
Appendix C. 
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Figure 3. Squirmer wake volume as a function of /3 = (4K,)/(3Va 2 ) where k is the dipole 
strength, V the velocity and a the radius of the squirmer. The wake volume is measured in units 
of the squirmer volume. 




Figure 4. Geometrical model of Rhodobacter sphaeroides, after Shum et al. (20101 



Having substituted the above calculations in the expression (3.9|, we readily obtain: 

v D = v (l/2 - /(/§)) . (4.2) 

The first term in the brackets is the added mass coefficient for a sphere Cd = 1/2. The 
function /(/3) is plotted on Figure [i] Thus, when |/3| < 1 the volume entrained by a 
squirmer equals that of a sphere of the same size moving with the same speed in inviscid 
fluid. For larger vd decreases and becomes negative for \f3\ > 2.4. 

4.2. A model for R. sphaeroides 

We now examine numerically the paths of passive tracer particles in the flow caused 
by a more realistic model swimmer, designed to represent the bacterium, Rhodobacter 



sphaeroides (Shum et al. 2010). The geometrical model consists of a prolate spheroidal 



cell body and helically-shaped flagellum of finite thickness, as depicted in Fig. [4] Each of 
these components is rigid and the flagellum rotates about its axis, which coincides with 
the major axis of the cell body. 



Employing the boundary element method described in Shum et al. (20101, we numer 



ically determine the flow field and swimming velocity of the swimmer which satisfies 
the Stokes equations, no-slip boundary conditions on the surface of the cell body and 
flagellum, and the constraints that zero net force and torque act on the swimmer. 
We begin with simulations of a bacterial swimmer with a cell body of length/ width 




Figure 5. (a) Deformation of an initially uniform rectangular grid of tracer particles as a 
bacterial swimmer approaches from infinitely far away, (b) Deformation (averaged over angles) 
of the material sheet. Plots of z-component of the net displacement as functions of initial radial 
position on logarithmic axes. Dotted and dashed lines are guides indicating the slopes of power 
law scaling. 



aspect ratio Ln/L± = 2 and a flagellum of helical pitch A = 2 and length L — 10. The 
non-dimensional lengthscales and timescales used throughout this section correspond 
respectively to lpm and 0.1 s for typical bacteria. Using the average flow field, we compute 
the trajectories of tracer particles initially arranged along the line segment z = 0, 0.1 ^ 
p ^ 30. The swimmer starts at z = —1000 and we run the simulations until it reaches 
z = 1000. The results are extrapolated in both directions to give the tracer displacements 
for an infinite swimmer path (see Appendix D for details) . We find that we can neglect the 
net radial displacement which is orders of magnitude smaller than typical displacements 
in the z-direction. 

Paths of tracer particles for different initial values of p are shown in Fig. [T] The distur- 
bance of the fluid can also be visualised as in Fig. [5ja), which illustrates the deformation 
of an initially uniform rectangular grid of tracer particles as a bacterial swimmer ap- 
proaches from infinitely far away. 

Double logarithmic plots of the net tracer displacements for an infinite swimmer path 
as functions of the initial tracer radial distance p are presented in Fig. J5jb). For tracers 
that start far from the swimmer the net displacement is negative and scales as p~ 3 in 
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Table 2. Comparison of Darwin drift for different swimmer geometries. Unless listed in the 
first column, parameter values used are: L\\/L± — 2, A = 2, L — 10, a = X/2n. 



Shape 



I'D 



I'D 



(from Eq. (3.7 1) (from simulations) 



Base 


-0.15 


-6.10 


-6.11 


L||/Ll =0.5 


-0.68 


-12.74 


-12.78 


L||/Lx =3.5 


-0.08 


-5.24 


-5.33 


L = 5 


-0.17 


-6.36 


-6.30 


L = 15 


-0.14 


-5.98 


-6.03 


A = 0.5 


-0.20 


-6.76 


-6.78 


A = 3.5 


-0.04 


-4.71 


-4.68 


A = 8,L = 20 


0.58 


3.04 


3.07 



agreement with Eq. (3.2). For those that start near to the swimmer the net displacement 
is positive, and diverges as p _1 as p — > 0. This is because a fluid particle ahead and 
precisely on the axis of the swimmer approaches the stagnation point at the leading pole 
of the cell body and is pushed indefinitely by the swimmer. 

We calculate the Darwin drift vd by a direct integration of the volume swept out by 
the material sheet for an infinite swimmer path, obtaining Vd = —6.29 in units of cubic 
swimmer body radii. Details of the numerical integration are given in Appendix D. To 



make contact with the expression (3.9 ), we first note that the swimmer had no discernable 
wake. Hence = v s = 4ir/3 s» 4.19. Q±, can be obtained by numerical integration of a 
moment of the force on the swimmer over the swimmer surface, 



1 



IzP 2 dS 



(4.3) 



giving Qj_/V = 0.15. Putting these numbers into Eq. (3.9), gives Vjj = —6.10, in good 
agreement with the value obtained by direct numerical integration. 

These calculations were repeated using the flow fields generated by swimmers with a 
range of shapes, but keeping the volume of the swimmer body constant. In particular, 
we started with the base parameter values, (Lt\/L± = 2, A = 2,L = 10), and varied each 
parameter in turn using the values: 

L\\/L± = 0.5,3.5, A = 0.5, 3.5, L = 5, 15. 

None of the swimmers had a wake and hence = 4.19 as before. The value of Q±/V 
for each of the geometries is listed in Table [2] The Table then compares the values of vn. 
calculated by using Q±/V in Eq (3.9) to that obtained by direct numerical integration 
of the volume of fluid displaced, showing pleasing agreement. The results indicate that 
for plausible swimmer shapes, the drift is of the order of the swimmer volume. It can be 



seen from Eq. (3.7) and Eq. (4.3) that the drift depends on the radial distributions of the 



propulsive forces on the flagellum and drag forces on the swimmer. The drift is enhanced 
by increasing the difference between the widths of the cell body and the flagellum helix. 
Swimmer shape is therefore critical in determining the Darwin drift. 
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5. Discussion 

The aim of our paper has been to describe the entrainmcnt of tracer particles by 
microswimmers. This is a step towards a full understanding of the ways in which bacteria 
can stir their environment, its likely importance in enhanced nutrient uptake and its 
possible relevance to oceanic mixing. Our results hold in the limit of zero Reynolds 
number, which is appropriate at the length scales of bacteria, and rely heavily on the 
observation that bacteria are force-free so that their far field is, to leading order, dipolar. 

The trajectories of tracer particles far from a swimmer that moves in a straight line 
along an infinite path form closed loops. This observation is at odds with experiments 
and simulations showing enhanced tracer diffusion in bacterial suspensions. Therefore 
we first present a full derivation of the tracer displacement which allows us to list the 
conditions necessary for closed loop tracer trajectories. We find that the loops are closed 
only if the tracer velocity can be neglected compared to the swimmer velocity. As the 
initial distance between the tracer and the path of the swimmer p decreases this condition 
is violated to an increasing degree, and a far field expansion shows that the net tracer 
displacement is negative and increases as p~ 3 for a dipolar swimmer, and as p~ 5 for a 
quadrupolar swimmer. 

Numerical results for a model swimmer show that as p decreases further the net tracer 
displacement changes sign. For tracers very close to the swimmer it becomes large and 
positive and diverges as p — > 0. This is expected as it reflects the behaviour of a tracer 
ahead and precisely on the axis of the swimmer that approaches the stagnation point at 
the leading pole of the cell body and is pushed indefinitely by the swimmer. However 
the total volume of fluid transferred by the passing swimmer, a quantity known as the 
Darwin drift, remains finite. We calculate the Darwin drift and show that it can be 



written simply as a difference of two terms, see Eq. (3.9 1. The first of these is, perhaps 
surprisingly, universal, as it depends only on the quadrupole moment and velocity of 
the swimmer. The second term has a pleasing physical interpretation, being equal to the 
sum of the volumes of the swimmer and of its wake. We show that the formula holds 
for squirmers, where the Darwin drift can be calculated analytically, and agrees with 
numerical results for a model swimmer resembling Rhodobacter sphaeroides. 

A second mechanism which leads to non-closed tracer trajectories, and hence tracer dif- 
fusion, is finite swimmer run lengths. Mixing enhancement due to random re-orientation 
of straight runs was first suggested by 



Lin et al. 


( 


2011) 


Polin et aL||2009 



flagellar desynchronisation that leads to a run-and-tumble manner of propulsion. It has 
been hypothesised ( Stocker & Durham 2009 ) that this feature has evolved as a stealthy 



strategy to avoid predators. The enhancement of biomixing may serve as an alternative 
evolutionary hypothesis. 

As a bacterial suspension becomes more concentrated, tracers will interact simultane- 
ously with several swimmers, and this will be a factor in determining their diffusivity. 
Moreover, in this article we have neglected Brownian motion. On average Brownian 
fluctuations do not destroy the loop-like topology of the tracer paths and a first approx- 
imation may be to add the Brownian and swimmer contributions to the tracer diffusion. 
However this remains to be verified. 

To summarise, we have given an account of the effect of entrainment on the paths of 
tracers moving near a low Reynolds number swimmer. This is one of the mechanisms, 
several of which are identified here, which may lead to enhanced tracer diffusivity in swim- 
mer suspensions. We hope that our results are a step towards determining the relative 
importance of various stirring mechanisms in different regimes of swimmer concentration 
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Appendix A. The scaling law A ~ p 3 

In this section we derive the scaling law for the passive tracer displacement as a single 
swimmer travels in a straight infinite path. We will assume the origin at the initial tracer 
location ry^ = and direct the z-axis along the swimmers's path, see Fig. [lj. 

When the particle displacement is much smaller than the swimmer displacement, 
|rx| <C |R|, the right hand side of the equation for the tracer motion ( |2.1[ ) can be 
expanded as: 



u(R(*)) - (r T • V) u(R(i)) + O ((|r T |/|R|) 2 ) . 



(Al) 



We solve the above equation using the perturbation theory. For the zeroth order term 

„(0). 



T ■ 



d.r 



(0) 



dt 



u(R(t)). 



(A 2) 



If the z-axis is directed along the swimmer's path, the equation (3.1| can be written 



as: 



u(R) = d z 



k\Js (R,k) 



In , Q±-Q\\ 



d z V s (R,k) 



0(R~ 



(A3) 



where R = |R| and d z — (k • V). Taking into account that 4f = Vd z we can integrate 



( A 2 1 along the swimmer location z and obtain: 



V 



U s (R(z),k) 



V R 



Q±-Q\\ 
2V 



9 2 U s (R(z),k) 



+ 0(R(z)- 3 ). 



(A 4) 



R(z) 



In determining the integration constant in the above expression we used the initial con- 
ditions Zi = — oo and r^ ^ = 0. It can be immediately seen that as the swimmer goes to 

infinity z +oo, |r^ | ~ R(z)^ 1 —> 0, i.e. the tracer trajectory is a closed loop, in agree- 
ment with the general result formulated in Section 2. Therefore the tracer displacement 
will be dominated by the first order perturbation Ar^" 1 . 
The equation for r^P reads: 



fir 



dt 



(4 0) -v) u(R(i)). 



(A 5) 



Substituting in the expressions for u and obtained above and keeping the terms up 
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to the order 0(i? -5 ), we obtain 



dr 



(i) 



dt 



V 

V 
k(Q± 



(U s • V) d z u s 



(u s • v) d z v 



On) 



2V 



3 z [(Us-v)d 2 Us] + o(ir 



(A6) 



The latter equation can be integrated along dz = Vdt. The term with the prefactor 
proportional to (Q± — Qy) is a derivative along the path and will, therefore, vanish upon 
integration from — oo to +oo. It is easy to see that the result has to be of the form: 

" 2 1 . _ «Qj 



1 



A-(i) _ r 

- ^V 2 p3 ' ~ Z p 4 



Op 



-5^ 



(A 7) 



where p is the distance between the swimmer's path and the initial position of the 
tracer, and Ci, C2 are (vector) constants. It turns out that the displacement in the 
direction normal to the swimming direction is negligible when the swimmer's path is 
infinite. Therefore we will focus on the tracer displacement along the swimming direction 
A = (k-Ar T ): 

A- C K " 1 +C KQ 1 

Here C\ and C2 are scalars. Let us calculate the numerical value of C\. According to 
(|A6|) and (TaI, 



Ct = 



/+00 
dz P 3 (u s -v)a,(k-u s ). 
-00 



(A 9) 



Detailed calculations show that C\ = n/16 ~ 0.2. We note in conclusion that when k = 0, 
the dominant contribution to A scales as p~ 5 . 



Appendix B. The Darwin drift 

B.l. Streamfunction 
For a cylindrically symmetric flow the streamfunction ip(p, z), defined as 

ip(p, z) = J —dp2npv z + dz2itpv p , (B 1) 

has a clear physical meaning: it equals the volume flux of incompressible fluid crossing 
a surface of revolution obtained by rotating about the z-axis a two dimensional curve in 
the coordinates (p,z). Without loss of generality we will assume that one of the curve 
ends lies on the z-axis and ip(p = 0, z) = 0. Incompressibility of the flow guarantees that 
ij) is a function of the coordinates (p, z) independent of particular curve choice. According 
to our definitions, for a uniform flow u p = 0, u z — —V and the streamfunction 

^ (p,z) = lvp 2 . (B2) 

In order to find the asymptotic shape of the streamfunction for a force-free swimmer, 
we note that in the frame co-moving with the swimmer with the speed V, v z = —V + u z . 



Then, according to (B 1 1, 

dij' 



. = -2npv z = 2irp(V - u z 
dp 
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Therefore, 



2n 



1 



dppu z 



(B3) 



where the prime denotes that the integral is taken along z = const. According to (3.1), 
the flow generated by the swimmer 



Q\\ 



d 2 



Q 2 v 2 | < 



o 



u Sz 



1 z 2 
r + ^ 

where r is the distance from the swimmer. Then is it straightforward to calculate that 



d z U Sz = r- 2 ((z/r) - 3(z/r) 3 ) , 
d 2 z U Sz = r- 3 (1 - 12(z/r) 2 + 15(z/r) 4 ) , 
V\U Sz = r~ 3 (1 + 6{z/r) 2 - 15(z/r) 4 ) . 
We obtain upon integration and a few straightforward transformations 



(B4) 
(B5) 
(B6) 



dppu z = K 



n 

p Z 



2 / 1- 3cos 2 i 

y ii ■ 



Q. 



1 + 3 cos 2 



\ -II 2 ' 2 

+ 0{r- 2 ) + C{z). 

Here C{z) is an arbitrary function at this point. In order to determine it we recall 

9-0 



(B7) 



dz 



lltpVp — 2-KpUp 



It can be verified that this expression holds if C'(z) = 0(z~ 4 ) for large z. Therefore 
C(z) = C\ +0(z~ 3 ). As far from the swimmer its influence on the flow must be negligible, 
we conclude that C± — and arrive at the final expression: 



2tt 



V 
~2 



KZ 



1-3 cos 2 



Q. 



l + 3cos 2 0\ 1 



+ 0(r- 2 ). 



B.2. Drift 



Following the original argument of Darwin (19531, we consider a circular material sheet 



of radius po placed in front of the swimmer at the initial distance Li, Figure [2j In the co- 
moving reference frame the incoming flow is uniform and has the speed V. If the swimmer 
did not perturb the flow, the material sheet would travel with the flow and after a time 
At assume the position z = Lf = Li — VAt. The fluid volume enclosed within the 
region swept by the material sheet, equals the total fluid volume through the boundary 
of this region during At, i.e. it equals ipo{po,Li)At, where ipo is the unperturbed flow 
streamfunction. 

When the incoming flow is perturbed by the swimmer, the total amount of fluid swept 
by the material sheet equals ip{pQ,Li)At. It can be seen from Figure [2] that 

ipAt + v D + v* -A = i/} At (B8) 
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where Vn is the total volume of ambient fluid entrained by the swimmer, is the fluid 
volume enclosed by a stagnation streamline around the swimmer, and A is the "excess" 
volume equal to 



A( Po ) 



dzw{p (z) - p ), 



(B9) 



where p(z) is the streamtube width. In the limit of large po and Li, tp — » and we 
obtain 



v D = A(p ) - u*. 



(BIO) 



Let us find the asymptotic shape of A(p ). We re- write the streamfunction (3.3) as: 



^ l+B(8) 3 -r(ff) 



where B(6>) = /3 sin 2 9 cos 9, T(9) = sin 3 9 1 



a 

-3 cos 2 



0((a/p)- 2 ), 



(Bll) 



, 1+3 cos-' 9 
2 ' '-L 2 



), /? = 2 K /(T/a 2 ), 



7|| = 2Q||/(Fa 3 ) and 7+ = 2Q_ L /(Va 3 ). Far from the body ip —> ip = nVp 2 , therefore 
on the streamline starting at p = po, z = 00, 



i>(p, z) _ 2 

Substituting ( B 12[ ) in (Bll) and denoting x — a/ po, y — p/po, results in 

y 2 (l + Bx 2 y- 2 - Tx 3 y- 3 + 0{p^)) = 1, p -»■ 00. 
After resolving this expression with respect to y we obtain 

y = 1 - ^B.t 2 + ^Tx 3 + 0(x 4 ), x^O. 



(B12) 



We are now ready to calculate the integral in (B9) in the limit of large pq. 



f 

A = np 3 I d(ycot9)(y 2 -l) 

i 

rf 



f d {l-\vx 2 + \Tx 3 + 0{x A ))cot9 



{-Bx 2 + Tx 3 + 0{x 4 )) 



f 



= Tzp% / d (cot 6») (-Bx 2 + Tx 3 + 0(x 4 )) 



irp 3 / d9 



P cos ftr — sin 9 I 7 



1 - 3 cos 2 9 1+3 cos 2 



+ 7-L' 



Til "t" T-L 

nf3a 2 po sin 6 1 + 7r J a 3 cos + n ; — — a 3 cos 3 9 



2 

-if 



x 3 + 0(x 4 ) 



2-KK 



Pq sin ( 



/ , 7T(Q||+Q ± ) 7T(Qx-Q||) 3 a 

cos 6* H — — cos 9 



V 



V 



+ OiPo 1 ) 



O^o 1 ). (B13) 



Appendix C. A squirmer wake 

We derive the expression for f(/3). For a squh 



mer 



^=p 2 



■£((?)'-(;)>•■ 



fl3 

r 
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On the stagnation streamline ip — 0. This occurs when 

r = a, (CI) 
i.e. on the squirmer surface, or on the surface 

-cos^ ri( :' + ri + 1) , r^r/a. (C 2) 

(30/2)( ri + i) 

This surface intersects the swimmer body at an angle 

cos6»* = 

Thus, the wake exists only if |/3| > 1. Let us calculate its volume. 

v wake = [ drd62Trr 2 smd = ^na 3 f0), (C 3) 

J wake " 

f(fi) = 2j dri^(l + coBd(n)), (C4) 



where 6(ri) satisfies ( C 2 1 and r*(/3) satisfies 

r^rl + r* + 1) - (3/3/2) (r* + 1) = 0. 



The integral in ( C 4 1 can be taken analytically: 



r 3 _ i 



31 , /r* + l\ rt ri r° 



log 



— r* 



30 ° V 2 7 2 3 5 

The dependence f(/3) is shown in Figure [ij 



Appendix D. R. sphaeroides simulation details 

Calculating the flow field: As we deal only with the case of a swimmer with an axisym- 
metric body in unbounded fluid, we need only compute the velocities once using the 
boundary element method. The velocities are constant in the reference frame of the flag- 



ellum and hence, we can calculate the mean direction and speed of propagation (Keller 



& Rubinow 1976). We choose the mean swimming direction to be along the z-axis and 
scale the spatial and temporal variables so that the volume of the cell body is 47r/3 cubic 
units, the volumetric radius is 1, and the mean swimming speed is 1. 

With typical geometrical parameters for the bacterial shape, we find that as the flagel- 
lum turns, the direction of the major axis of the cell precesses with a negligible amplitude 
about the direction of mean swimming. We therefore make the approximation that the 
swimmer is orientated in the z-direction at all times. We find that the period of the 
turning flagellum is relatively short; more than 10 complete turns are required for the 
swimmer to advance by one body length. We therefore consider the motion averaged over 
a period of revolution. 

We tabulate the flow velocity on a grid of ~ 35,000 points. The grid spans the ranges 
— 1000 ^ z ^ 1000 and ^ p ^ 150, and a higher density of points is used near the 
swimmer to more accurately capture the variations in velocity. Between grid points, we 
use a cubic spline interpolation to approximate the average flow field. 

To characterise the effect of a bacterial swimmer on the surrounding medium, we track 
a number of points in the fluid that move with the local average flow field, using the 



MATLAB ODE solver, odell3 to integrate Eq. (2.1 1 in time. 



18 



D. O. Pushkin, H. Shum and J. M. Yeomans 



Extrapolation to an infinite swimmer path: We begin simulations at time t = —1000, 
which means the swimmer is at position z = —1000 and the trajectory of each tracer 
particle is computed until the swimmer passes the particle and reaches the distance 
z = 1000. To extrapolate the net tracer displacements to an infinite swimmer path we 
note that, for a dipolar far field flow, the tracer positions will follow the limiting behaviour 

U-cc) + C*_Jt, t->-oo, 
W \z{oo) + C^/t, t-too. 

A nonlinear least squares algorithm is then used to estimate the limiting positions z(— oo) 
and 2(00), and hence the net displacement along the z-direction, A = 2(00) — z{— 00), 
for an infinite swimmer path from the finite computed trajectories. 
Calculating the entrained volume: Net tracer displacements are obtained numerically for 
initial radial tracer positions 0.1 < p < 30. The curve is extrapolated to all values of p 
using the ansatz 

Mp) = 1 C0/P / \ 0< ' <(U ' (D2) 

(see Fig.^b)). The curve of net tracer displacements A can be divided into a negative 
region and a positive region with the transition at the point p = p , say. Obtaining p by 
interpolation of the numerical data, we calculate the volumes of the bodies of revolution 
in the two regions of positive and negative net tracer displacements by evaluating the 
formulas 

rpo roo 

v = 2tt P A(p) dp, Voo = / 27rpA(p) dp, (D 3) 

JO J on 



using quadrature within the range 0.1 < p < 30 and analytical integration of (D2) for 
the contributions outside this range. The numerical results show that typically, a volume 
of fluid about twice that of the swimmer is pulled forward with the swimmer across any 
given interface perpendicular to the swimming direction. At the same time, fluid further 
away from the path of the swimmer is pushed backwards. The net flux is the Darwin 
drift, given for different swimmer shapes in Table [2] 
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